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When two types of side chains (A,B) are densely grafted to a (stiff) backbone and the resulting 
bottle-brush polymer is in a solution under poor solvent conditions, an incompatibility between 
A and B leads to microphase separation in the resulting cylindrical brush. The possible types of 
ordering are reminiscent of the ordering of block copolymers in cylindrical confinement. Starting 
from this analogy, Leibler's theory of microphase separation in block copolymer melts is generalized 
to derive a description of the system in the weak segregation limit. 

Also molecular dynamics simulation results of a corresponding coarse-grained bead-spring model 
are presented. Using side chain lengths up to N = 50 effective monomers, the ratio of the Lennard- 
Jones energy parameter between unlike monomers (eab) and monomers of the same kind {eAA ~ 
^bb) is varied. Various correlation functions are analyzed to study the conditions when (local) Janus 
cylinder-type ordering and when (local) microphase separation in the direction along the cylinder 
axis occurs. Both the analytical theory and the simulations give evidence for short range order due 
to a tendency towards microphase separation in the axial direction, with a wavelength proportional 
to the side chain gyration radius, irrespective of temperature and grafting density, for a wide range 
of these parameters. 



1. INTRODUCTION 

Enabled by progress in chemical synthesis (see [l|, 0| 
for reviews), macromolecules with "bottle-brush" archi- 
tecture where flexible side chains are densely grafted to 
a backbone, have found much recent interest (e.g. [ij- 
[isjV Such molecules may be useful for various appli- 
cations (such as sensors, actuators, building blocks in 
supramolecular structures 0^1l), since these systems are 
stimuli-responsive polymers, exhibiting large conforma- 
tional changes when external conditions vary. Apart 
from these synthetic bottle-brush polymers, also biopoly- 
mers with a related architecture are abundant in nature, 
e.g. proteoglycans 1J|. These brush-like polymers con- 
tain a protein backbone with carbohydrate side chains, 
and are held responsible for a large variety of biologi- 
cal functions (cell signaling, cell surface protection, joint 
lubrication, etc. |15l4l7|). At the same time, the in- 
terplay between entropic effects and various enthalpic 
forces in these "soft" objects makes the understanding 
of structure-property relationship of such bottle-brush 
polymers a challenging problem of statistical thermody- 
namics [18.-55.]. 

In the present paper, we focus on the interplay be- 
tween chain conformations and local order in binary A,B 
bottle-brush polymers, from the point of view of both 
theoretical arguments and Molecular Dynamics Simula- 
tions. Previous work [s^ 41 1 has focused on the pos- 



sibility of microphase separation in the form of "Janus 
cylinders" i.e., the cylinder splits in two halves (in 
the most symmetric case, where the lengths Na,Nb of 
both types of side chains and their grafting densities 
aA,<yB are equal), such that the A-B-interface contains 
the cylinder axis (taken to be along the z-axis hence- 
forth). Then it was pointed out I46| that due to the 
quasi-one-dimensional character of this ordering, for fi- 
nite cylinder radius R, no true long-range order of Janus- 
cylinder-type should be expected: rather there exists a 
finite correlation length over which the orientation of the 
interface plane decorrelates. In addition, it was specu- 
lated 47] that also the ratio eab/^ between the interac- 
tion strength eAB of unlike monomers and monomers of 
the same kind (eaa — ^bb — ^) matters: if the formation 
of A-B interfaces is energetically much more unfavourable 
rather than the formation of polymer-solvent interfaces, 
a phase separation along the axial (z) direction of the 
cylinder could occur also in the form of a double cylinder 
(the A-rich polymers form a separate cylinder from the 
B-rich ones, both cylinders touch each other along the 
z-axis) [47]. Also intermediate cases ("Janus dumbbell" - 
like cross sections of the bottle-brush) were suggested 14711 
and some indication of such structures were observed [55| . 
However, again only finite correlation lengths of such or- 
derings along the z-axis can be expected. In addition, 
one needs to consider that for not very large grafting 
densities the bottle-brush polymer under poor solvent 
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conditions does not form a structure that is homoge- 
neous along the z-axis (the direction of the backbone; 
note that we restrict attention to rigid backbones only, 
for simplicity) , but rather inhomogeneities in the form of 
a "pearl- necklace" -structure form 43, 51, 56| . 

The problem which we investigate in the following sec- 
tions is whether it can be favorable to form a microphase- 
separated structure that is inhomogeneous along the z- 
axis. Following the information from experiment that one 
side-chain can be grafted per backbone monomer [lo| , we 
assume an alternating grafting of the two types (A and B) 
of side chains [55[- Of course, this situation is essentially 
equivalent to the case where symmetrical block copoly- 
mers 57h61| would be anchored with their A-B junction 



to the backbone. It is then also interesting to ask how 
such a situation compares to the case where such block 
copolymers fill a cylinder (at the same density) with- 
out constraint on the locations of the A-B junctions. Of 
course, the latter problem has already been studied in 
other contexts 62h76I| . Here we generalize the treatment 
of Ref. [iHl to take this constraint on the location of the 
A-B junctions into account to formulate a weak segrega- 
tion theory 77], |78| of ordering in binary bottle-brushes 
(Sec. 2). Then we present data from Molecular Dynamics 
simulations, extending the study of Ref. (ssj l by consider- 
ing now the variation of properties with eAs/e (Sec. 3). 
We summarize the main findings of our work in the con- 
cluding section (Sec. 4). 



2. WEAK SEGREGATION THEORY OF 
MICROPHASE SEPARATION IN 
TWO-COMPONENT BOTTLE-BRUSHES. 

2.1. THE CONDENSED STATE OF THE 
TWO-COMPONENT BOTTLEBRUSH. 

The general idea of our theoretical approach is as fol- 
lows. We consider the case when incompatibility of both 
sorts of blocks with the solvent is high enough to form 
a condensed state (globule) , which in the case of a rigid 
backbone would acquire the form of a cylinder uniformly 
filled by the monomers of the side chains. The height 
H of the cylinder would be just the length of the rigid 
backbone and the radius R would be related to the equi- 
librium monomer density po inside the cylinder. Such 
a condensed state could be considered as an equivalent 
of a bottlebrush compressed up to the monomer density 
po via squeezing it into a tube of the radius R. To con- 
sider the composition fluctuations in such a condensed 
(squeezed) bottlebrush we employ some new ideas, which 
were put forward in refs Tl], 75 1 to generalize the semi- 
nal Leibler theory 77 1 (see also refs 57 -gH) to confined 
polymer systems. To begin with, we find, based on the 
theory of polymer globules by I. Lifshitz et al. [t^, the 
value of po, which is to provide the minimum of the two- 



component bottlebrush free energy given the values of 
the chi-parameters Xij (* similar approach has been em- 
ployed by Borisov et al. [i^ for one-component bottle- 
brushes in a poor solvent). To determine the equilib- 
rium value of a two-component bottlebrush forming a 
condensed state due to its overall incompatibility with 
the solvent we notice that, according to I. Lifshitz et al 
[l^, the free energy of any fiexible polymer system is 
(in a zeroth approximation) just the sum of a structural 
(entropic) contribution, which allows for the effect of con- 
nectivity (" linear memory" in Lifshitz terms) and an en- 
ergetic contribution, which allows for the van-der-Waals 
attraction and the excluded volume repulsion effects: 



F = F3,,({p.(r)})+F* ({p,(r)}) 



(1) 



The precise form of the structural free energy 
Fstr ({Pi(r)}) as a functional of the spatially non-unform 
distributions {pi(r)} of the local number densities of the 
repeat units of the i-th sort (z = B) is determined by 
the actual microscopic structure of the bottlebrush. We 
analyze it somewhat later when studying the bottlebrush 
stability with respect to its longitudinal, angular or he- 
lical ordering. But now, when we focus on the basically 
uniform cylindrical globule, the structural term could be 
simply disregarded as compared to the energetic contri- 
bution to the free energy: 



F^F*{{p,{v)})^TVr {{-Pi\) 



(2) 



Here V — ttR^H — 2MTTR'^a is the volume of the cylin- 
drical globule (where 2M and a are the total number of 
both A and B side chains and the distance between the 
neighboring A and B chains), Pi is the value of the lo- 
cal number density of the repeat units of the i-th sort, 
averaged over the whole volume V of the bottlebrush, 
and the specific (per unit volume) free energy /* ({p^}) 
depends on the type of the interactions between the sol- 
vent molecules and repeat units (monomers) forming the 
bottlebrush. For simplicity, we assume in this paper that 
these interactions correspond to the conventional incom- 
pressible Flory-Huggins lattice model, in which case 

vof *{<!>) = (1-0) In (1 -(/.)+(/. 

+ XAB(I)A(I)B +XAs4>A(t>S + XBS(t>B(l)S, (3) 

where vq is the excluded volume assumed, for simplicity, 
to be the same for both sorts of the repeat units and 
solvent molecules, (r) = vopi (r) is the local volume 
fraction of the particles of the i-th sort within the globu- 
lar bottlebrush, (j) — (j)^ + (t>B is the total polymer volume 
fraction (volume fraction of all monomers), (/)5 = 1 — is 
the volume fraction of the solvent molecules and Xij ^^e 
the conventional Flory-Huggins energetic parameters. 

For the homogeneous bottle brush we are considering 
here, (jji (r) = (pi = fi<j) = i^oPi, where fi is the fraction 
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of the monomers of the i-th sort, and, thus, the function 
/* {{Pi}) can be rewritten in the form 



vof* ({Pj) 



vof* (0) 
(1-0) ln(l 



+ (j)-xr + x(^ (4) 



where x = XAsfA + XBsfs-XABfAfB is the effective X" 
parameter describing the overall bottlebrush-solvent in- 
compatibility and the parameter x = xasIa + Xbs/b 
does not affect the equilibrium value of (f). 

Indeed, taking into account that the polymer volume 
fraction within the uniform bottlebrush is 



(5) 



where N ~ Na + Nb is the total average number of all 
monomers A and B per the distance 2a of the backbone, 
one can reduce the total free energy ([2|) to the form 



F = TMN 



)ln(l 



+ - X 0^ + X ( 



(6) 

Minimization of the free energy ([6]) with respect to 
results, finally, in the desired equation (it corresponds 
to the so-called volume approximation of the polymer 
globule theory [7§|) for the equilibrium polymer volume 
fraction of the two-component bottlebrush: 



voP*/T=-h-i{l-cj))~cj)~xc^^ 



0. 



(7) 



Obviously, the finite equilibrium polymer volume frac- 
tion of the two-component bottlebrush which is the 
solution of eq Q, depends on the value of x only and 
has a physically meaningful positive value only for 



r = X- 1/2 > 0. 
The corresponding asymptotics read 



= 3t - (27/4)t2 + 
= 1 - exp (1 + x) , 



T < 1 

X> 1 



(8) 



(9) 



It is worth to notice that in so-called Hildebrand ap- 
proximation [soj . in which the value of the x-parameter 
for the i-th and j-th species is proportional to the squared 
difference of their solubility parameters: 



Xij = vo (St - Sjf^ (2r) , 



the effective x-parameter is always positive: 



X 



{2f-l + xf 



-XAB, 



where x = (26s — Sa — Sb) / {Sb — Sa) is the selectivity 
parameter |8l| . If the A and B monomers are compatible 
{6a = 6b) then x = XAs = Xab- 

To conclude this subsection we discuss the conditions 
of validity of our consideration. First, it is worth to notice 



that the bottlebrushes under consideration are quasi-lD 
systems. Therefore, their equilibrium volume fraction 
stays finite even in the 0-solvent, where ^ NP (here 
I is the statistical segment supposed, for simplicity, to be 
the same for both sorts of the side chains). Indeed, it 
follows from eq ([5]) that in this case 



(f)e ^ Li/ (tto) , 



(10) 



where we introduced the reduced distance a = A/l be- 
tween the neighboring side chains and the Lifshitz num- 
ber Li = vo/l^, which characterizes the chains' flexibility. 
Comparing eq and the first of eqs ^ and taking into 
account that the condensed bottlebrush should be much 
denser than that in 0-solvent we conclude that the first 
of the validity conditions reads 



K = 0e/0 = R^/{NP) w Li/(37raT) < 1 , 



(11) 



On the other hand, the volume approximation (jT]) is 
known [79] to be valid when 



i? > Tc 



(12) 



where Tc is the correlation radius within the globule. For 
dense globules (1 — (^^lorr^l) the condition holds 
always but for semidilute globule (/3/3 ^ t <C 1) the cor- 
relation radius is known 7£ 
of r: 



831 to increase with decrease 



l/r. 



(13) 



It follows from eqs ([5]), ([9]) and (fT3|) that close to the 
0-temperature the condition (|11[) is satisfied only when 
the side chains are long enough: 



iV» l/((/)eT). 



(14) 



Condition is rather severe since both quantities (pQ 
and T are small. Thus, the fact that the monomer density 
profile usually observed in computer simulation of bot- 
tlebrushes in poor solvents close to the 0-temperature is 
not step-like, which is expected for a condensed (glob- 
ular) state, but rather smooth [Hll, [H^ is explained by 
insufficiently high degrees of polymerization of the side 
chains, which do not satisfy condition (fT4|). 



2.2. THE RANDOM PHASE APPROXIMATION 
FOR CONFINED POLYMER SYSTEMS. 

General theory. 

To consider the fiuctuations within the two-component 
condensed bottlebrush we are to expand the total free 
energy ([1]) in powers of the density fiuctuations counted 
from the uniform globular state (5'i(r) — pi{v) — 'Pi)'- 



F = Fo + AF2 + AFg -I- AF4 



(15) 
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where Fq is defined by eqs ([2]) and ([6]), the quadratic term 
of the expansion reads 



AF2 = |y" r,,(ri,r2)*,(ri)*,(r2)dVidl/2 



(16) 



(in eq (jl6l) and thereafter we employ the rule of summa- 
tion over the repeating indices) and the next terms of the 
expansion (fT5|) are defined similarly fll\ . 

Our purpose in this subsection is to find the explicit 
form of the kernel r(ri, r2) and conditions ensuring that 
the quadratic form AF2 is positive definite, which implies 
that the fiuctuations are small. It follows from eqs ((Ij - 
and (HH) that 



Tapir I, r2) 



T-^5^F{{p,ir)}) 



= 7Q/}(ri,r2) - CQ/}(ri,r2). (17) 
Here we introduced new designations 

({P.(r)}) 



7a^(ri,r2) 



(18) 



Ca/3(ri,r2) 



T-^S^F* ({p,(r)}) 



Spa{ri) 5pi3{V2.) 

where the number matrix C — ||Cq^|| reads 



5{yi - T2) Cap, 
(19) 









dpa dpp 



2XAS ■ 

k - 6-^^ 



k-(t). 



2x 



BS 



(20) 

and k = XAS + Xbs - Xab- 

Thus, our problem is reduced to calculating the second 
derivative of the structural free energy functional. 
For this purpose we consider an auxiliary thermodynamic 
potential 



$str iW^ (r) , T}) - -T In Z ({^, (r) , T}) 



(21) 



which has meaning of the free energy of an ideal polymer 
system with a specified architecture (in our case it is 
the two-component bottlebrush) affected by a set ipi{r) 
of external fields applied to the particles of the i-th sort. 
The thermodynamic potential $str {{fi (r) , T}) could be 
readily calculated (see below). On the other hand, it is 
directly related to the desired thermodynamic potential 
_F'str(pj(r)). Indeed, the partition function Z appearing 
in eq (|2ip is the integral over all possible non-uniform 
density distributions pi (r) of the particles of the z-th sort 
in the volume of the system: 

^ = -J Wp^i^) (- i^str ({p.(r)} , (r)}) / T 



(22) 



with 



^^str {{p^{r)} , W^ (r)}) = F,,, {p,(r)} + / p,(r) ^,(r) dV. 



Definition and calculation of the density integral Z is, 
generally, rather cumbersome but it becomes trivial in 
the saddle-point approximation: 

$str iW^iv)}) = min F,,, {{p,{v)} , {ifi, (r)}) (23) 



It follows from ((23|) that if the thermodynamic poten- 
tial -Fstr ({pj(r)}) is known then the thermodynamic po- 
tential $str iWii^)}) is parametrically defined as follows: 

ipair) = - SF,^, i{p^{r)})/ Spair) (24) 
'i'str iWdr)}) ^ F,,, ({P.(r)}) 

SFstr {{p^{r)}) , , . 

P^'^^) — X — 7:r\ — (25 

dpa[r) 

Eqs (|24|) and ([25]) imply that the thermodynamic poten- 
tial $, which is a functional of all external fields (pi{r) 
applied to the particles of the system, is the Legendre 
transform of the thermodynamic potential F, which is a 
functional of all number densities Pi{r) of these particles. 
Therefore, the reciprocal relationships hold: 



Pair) ^ S^striWtir)})/ Sipair) 

Fstr {{p^{r)})^^s^r{W^ir)}) 

J 0(pa{r) 
It foUows from eqs (HH) and ([241) -([21]) that 

6'F^,A{P^ir)}) 



7a/3(ri,I"2) = 



Spa{ri) Spfj{r2) 
5ipa{ri) 



(26) 



(27) 



(28) 



[r) — const 



pj (r)— const 



'5p/3(l-2) 

Furthermore, if we define the structural matrix 

<^2$str({^.(r)}) 



ffa/j(ri,r2) 



Sipairi)Sipi3{r2) 
Spa{ri) 



(29) 



V(r)=0 



V(r)=0 



6(pf3{{r2)) 
we obtain 

lapiri^r) gp^{r,r2)dr = Sa-ySiri - rs) 



(30) 



Calculation of the structural matrix for copolymers 
in bulk. 



The derivation of the explicit expression for the struc- 
tural matrix g of copolymer systems under confinement 
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is a straightforward extension of that in bulk. So, we re- 
mind the reader how the structural matrix is derived in 
bulk. 

Instead of calculating the density integral ([22l) we con- 
sider an equivalent expression for the partition function 
Z{{(fi (r)}) of an ideal polydisperse n-component system 
of linear macromolecules affected by the external fields 
{(fii (r)} applied to the monomers of the i-th sort. The 
partition function Z{{ipi (r)}) can be written directly in 
terms of the corresponding discrete microscopic model: 



(31) 



Here 5* is the macromolecular structure, which for linear 
chains is just a sequence of integers a (i) equal to the 
number a of the sort of the z-th monomer, Ms is the 
number of chains with structure S in the whole volume V 
of the system and Zs {{(pi (r)}) is the partition function 
of the chain S: 

Zs (r)}) = JdTsF {Ts) exp (^'jy^i^) (r.) /T^ 

(32) 

where Ts = (ro,..,rAr) is a point of the configuration 
space of the chain S and 



N-l 



i=0 

is the distribution function of the ideal chain S, the ef- 
fect of connectivity being taken into account by the bond 
functions g (r) , which describe the correlations between 
the neighboring monomers of the chains. For simplic- 
ity, we assume that the bond function does not depend 
on the sorts of the neighboring monomers and has the 
conventional normalized form 



g(r) = exp(- ir^ (21^)) / {27rf / 3) 



3/2 



(33) 



where I is the statistical segment length. Thus, in the 
thermodynamic limit 

hm Zs{WAv)})l^^o/y='^- 
The thermodynamic potential $str ({'/^i(j^)}) reads 

$str ({^.(r)}) = 5]M5<i>s (34) 



structure implies that the matrix ([29)1 is a properly 
weighted sum of molecular response functions: 

5a/3(ri,r2) = ^ps5i^^(ri,r2). 



Here ps = Ns/V is the average number density of the 
macromolecules S and 



S'Zs iW.ir)}) 



Vi(r)=0 



where we employed the definitions (|32p and ([55]) . Sum- 
ming the second field derivatives of Zs {{(p{r)}) and col- 
lecting similar terms one gets 

5i5(ri,r2) = m(f) S (n - r,)+Y,,.%^ (n) g(")(ri,r2) 

n>l 

(36) 



where 



(ri,r2) = J ff(ri-r)g(")(r,r2)dr, (37) 



mi^'' is the number of the monomers of the i-th sort, 
which belong to the macromolecules S, and v^p{n) are 
the numbers of pairs of monomers of the a-th and /3-th 
sorts, which are separated by n bonds. The structural 
matrix, which was first introduced in the form p6p by 
one of us [^^,^] (see also 85] and references therein), 
is a discrete analog of the expressions for the correlators 
obtained in the continuous limit (like the Debye function 
for homopolymer and the Leibler expressions [t?! for di- 
block copolymers). 

To calculate the structural matrix explicitly we rewrite 
it as the inverse Fourier transform: 



J exp {iq (n - r2)) gg^ (?) 



(38) 



where the direct Fourier transform 



9a/3il)^V ^ J 5i'^^(i"i,r2)exp(iq(r2 -ri)) dridr2 



^i'^S^p + Y,,.%\n) g-{q) 



(39) 



TI>1 



is an algebraic function of the Fourier transform of the 
bond function: 



where 



$s [W^ {v)}) = -T\nZs{W, (r)}). 



(35) 



Additivity of $str with respect to the terms 

characterizing macromolecules with different chemical 



5 ((?) = dVg (r) exp (zqr) = exp {~Pq^/&) . (40) 



The explicit finite expressions for homo- and some 
block copolymers can be found in refs 83, 84j and 
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Calculation of the structural matrix for confined 
copolymers. 

Now, we make use of the fact that and exp(iqr) 
are the eigen values and eigen functions of the Laplace 
operator —A acting in the infinite space. The natural 
idea 71, 7^ to generalize expression (l38t for the struc- 
tural matrix to confined copolymers is to rewrite it in the 
form of a somewhat more general series 

5aMri,r2) = Y^'^ (^i) Y^^) (41) 

s 

where now {A^} and ^jf^ (r)| are the full sets of the 

eigen values and the orthonormalized eigen functions of 
the Laplace operator —A acting in a finite volume un- 
der specified boundary conditions and g (A) is a matrix- 
function of the parameter A. 

Then the kernels 7a/3(ri,r2) and Cq,^ (ri,r2) read 

7a/3(ri, r^) = (g"' ^^'^ (ri) ^^'^ (^2) , (42) 

S 

ca^i-i, r2) = C^pY. ^^'^ (^i) ^^'^ (^2) , (43) 

s 

which results in the following final expression for the de- 
sired kernel rQ^(ri, r2): 

r.Mri,r2) = Y ^^'^ (^1) (G"' i>^s))^,Y(^^ (r^) (44) 

S 

with 

G-i(A)-g-i(A)-C. (45) 

which is a generalization of the well known RPA equation 
for polymers in bulk: 



(46) 



We conclude that the condition ensuring that the 
quadratic form Ai^2 is positive definite is that ensuring 
that the matrix G has positive eigen values for any al- 
lowable value of A. 

In particular, for confinement in a cylindrical capillary 
of the radius R the eigen values and the eigen functions 
of the Laplace operator read: 

-X = q' + (C/i?)', (r) = J„ ^) e^^'^-""^ , 

(47) 

where r =(a, r, z) is the radius- vector in cylindrical coor- 
dinates, A/"^ — ^ /q^ i^riTi) ^ dr^ is a normalization 
factor and s™ is location of the m-th cxtremum of the 
Bessel function J„(a;) of the order n. 

Accordingly, the quadratic form (|16p takes the form 



T 



n=Or?j.= l •' 



dq 



where the coefficient ^^(g, n,m) reads 



For future comparison with ordering in two-component 
bottlcbrushes under investigation we present here the ex- 
plicit form of the structural matrix for diblock copoly- 
mers AnBjn- 



9'r^''is^,q)^cj,Nfo{Q'x^) 



(49) 



free 

9ab 



{s';^,q) = cl,N^{Q'xA)ij{Q^XB) 



Here N ^ n + m, xa = I - xb = n/N, Q'^ ^ k"!^ + Q'^, 
Q2 = q^pN/6, C = K'l/Rf N/6 = (s™)' /(6k), C is 
location of the m-th extremum of the Bessel function 
Jn{x) of the order n, parameter k is defined by eq (jlip 
and we introduced functions 

r I ^ exp(^j/) ~l + y (1 - exp(-y))2 

jD[y) = 2 2 ' ^{y) = 2 ■ 

(50) 



2.3. THE RANDOM PHASE APPROXIMATION 
FOR TWO-COMPONENT BOTTLEBRUSHES 

The derivation presented above for the confined block 
copolymers can be now extended to the tethered copoly- 
mers. For this purpose it is sufficient to replace the ex- 
pression (I3ip for the partition function of the correspond- 
ing confined system by that of the tethered system under 
consideration. In this paper we consider for definiteness 
the two-component bottlebrush consisting of 2M pairs of 
A and B side chains regularly distributed along the back- 
bone. The partition function of such an ideal bottlebrush 
affected by the external fields {(pi (r)} , reads: 

M-l 



II Zn^ {{ipA (r)} , 2na) x 

71=0 

M-1 

n^iVB({V'B(r)},(2n + l)a) (51) 



where M is the number of side chains of each sort, 
Zn {{^ (r) , zo}) is the partition function of an A^-mer, 
which is attached by one of its ends to the rigid backbone 
aligned along the z-axes at the point zq and affected by 
an external field (p (r) : 



N 



Zn (i")} : ^o) = n ( '^^^ 



y^(r») 

T 



g (r,_i - r, 



X (G i(g,n,m)) 



(48) 



(52) 

Here ro — (0, 0, zq) and g (r) is the bond function defined 
above by eq ([55)1 . the statistical segment length I being 
the same for A and B chains. 
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Then the thermodynamic potential ^str {{fii'r)}) 
reads 



M-l 



n=0 
M-1 

+ Y.'^^B{{VB{r)},r^) (53) 



where = (0, 0, 2na), r,f = (0, 0, (2n + 1) a) and 

$jv iW (r)} , zo) - -Thi Zn iW (r)} , zq) • (54) 

Additivity of $str ({v'i(r)}) with respect to the terms 
depending on different external fields implies that the 
matrix (j29p of the molecular response functions contains 
diagonal elements only: 



gAA{rur2)= ^ 



n=0 

M-1 .2 



Sip{ri) dip{r2) 



'^' S'^^, ({^(r)},0 



n=0 



V(r)=0 



V(r)=0 



5AB(ri,r2) = 5BA(ri,r2) = 0. 
On the other hand, 

<52$jv(Mr)},^o) _ <52^Ar({^(r)},zo) 



(55) 



(5(y3(ri)(5.^(r2) 



^(r)=o <5<p(ri)(5^(r2) 



V(r)=0 



<5Zjv({¥'(r)},zo) 5ZAr({<^(r)},zo) 



(5^(ri) 



(5(/3(r2) 



(56) 



¥'(r)=0 



where we employed the definitions ([5^ and ([5i)) . Sum- 
ming the second field derivatives of Zjv ({(p(r)}) and col- 
lecting the similar terms one gets 



gn{vx,V2) = Ai (ri)(5(ri - r2) + Bi (ri,r2) 
- d (ri,r2) , 



where 



M-l AT 



(57) 



(58) 



n=0 ; = i 



Jlf-1 AT 



i3,(ri,r2) = ^ 5^5^') (rj,,ri) g^") (ri,r2) 

n=0 i = l 

+ EE3^'Ur™,r2)5^™)(r2,ri) (59) 



A/-1 JV 



AT 



C,(ri,r2) = EE^^'^ ('•-'•i) Eff^"^ (^"''•2) (60) 



and the kernel g^^^ (ri,r2) is defined by eq ((37)) . 

Based on the previous consideration, we should use for 
^(fe) (ri,r2) the representation 



00 OC /"J 

5^'^^(ri,r2) = EE/|-P 

n— m— 1 



/cA (q, n, to) 
6 



xy,r(ri)y-(r2) 



(61) 



with the eigen values A and eigenfunctions F,™ (ri) de- 
fined by eq (|47l) . 

However, unlike the situation with free (untethered) 
copolymers under confinement, the correlators (|57p for 
bottlebrushes are, in general, not reducible to the form 
(|4Tt . Accordingly, the quadratic form (|T6)) does not take 
the desired diagonal form (j48|) . 

Fortunately, it could be shown that for the condensed 
state the correlators can be expanded in powers of the 
small parameter k defined by eq (jlip . Therewith, keeping 
only the dominant (zeroth approximation in k) term of 
this expansion one arrives precisely at the diagonal form 



^■^2 = -^E E / 2^^"(^' ^'^(^' 

n=Or?x=l 

X (G"^(g,n,m))^^, 

with 

G"""^ ((7, n, to) = g"""^ ((7, n, to) — C , (62) 
where the matrix C is defined in eq (I20p and 

5ao(«,C9)=0A^/i?(Q'2;„) (63) 



n = 



gt^n,s^,q)^<i>N jnyO'x^), n>\ (64) 
5ab - 0, n > 0, (65) 

with all parameters and functions defined after eqs (|49p . 



2.4. CORRELATIONS AND ONSET OF 
ORDERING IN TWO-COMPONENT 
BOTTLEBRUSHES: THEORETICAL 
PREDICTIONS. 

The general condition of the multi-component block 
copolymer disordered state stability with respect to mi- 
crophase separation is well known 78] . For block copoly- 
mers confined in a cylindrical tube and the condensed 
two-component bottlebrushes it should be reformulated 
as follows: 



n=0 i=\ 



m—l 



A = min A.; (n, to, q) > 0. 



(66) 
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where (n, m, q) are the eigen values of the correlation 
matrix defined by eq ((62)) given the values of the 
"quantum numbers" n and m characterizing the types 
of the angular and radial concentration modulation and 
wave number q characterizing the wave length of the lon- 
gitudinal concentration oscillation (75| , and the minimum 
is to be sought on the set of all possible values of n and 
771 and on the half-axes q >0. 

Accordingly, the condition, which determines the sur- 
face (in the space of the structure and energetic parame- 
ters of the system under study) where the uniform state 
of the latter loses its stability (spinodal) , reads 



. The corresponding eigen values are 



A{n*,m*,q*) = 0, 



(67) 



where we label with asterisks those values of n,m,q that 
actually satisfy eq (1571) . 

Based on the WST experience [stI - IHH one can arrive 
readily at the following conclusions, i) The value of n* 
defines the type of the angular modulation of the mor- 
phology arising after the uniform state looses its stability. 
In particular, if 7i* = then a radially symmetric mor- 
phology arises, if tt,* = 1 then the Janus-cylinder-type 
morphology would occur, \i n* — 2 then the cross-like 
ABAB morphology would arise etc. ii) If g* = then 
the arising morphology is expected to be homogeneous 
along the backbone direction, otherwise a ID morphology 
periodic along this direction would occur with the period 
L « 'iiT/q* . 3) the value oim (more precisely, s^) defines 
the type of the radial modulation of the arising morphol- 
ogy. (It is worth to notice that such a straightforward 
extension of the results of the bulk stability analysis to 
the quasi-one-dimensional ordering in tubes and bottle- 
brushes is only approximate [ig, 75| as mentioned above. 
We address this issue in more detail elsewhere.) 

To demonstrate the general results let us compare or- 
dering in bottlebrushes and free diblock copolymers con- 
fined in a cylindrical tube. It follows from eqs (PH)) and 
(l62l) - dMl) that for bottlebrushes 



{9%y^-2xBS 



(68) 



whereas for free diblocks in a tube 

/ {^'-''Taa-'^Xas (g 
(^oG)-i = 1^ 



' AB 



with the components of the matrix g^"^"" defined in eq 
(SSI), k = XAS + XBS - XAB and XaS = XaS - (20s)"^ 
For symmetric blocks the eigen vectors of the matrices 
(pl and dMl) are e_ = (l,-l)/v^ and e+ = (1,1)/V^ 
that correspond to the order parameters 

*_(r) - {S4>A {v)-5^B 

^+ (r) = {64>A (r) + 8ct>B (r)) 7^2 = -54>s (r) /V2 



h}^y (77, 777, q) = (5^b^ in, g)) - XAB. (70) 
K'^^ (77, m, q) = {g'fj, (n, s™, q))'' - X+, (71) 



for bottlebrushes (with x+ = ^XAS ~ XAB — '^4>s^) ^^^^ 

A^I- (s™, q) = {g'^A {^n,q) " 9ab (C, 9)) - XAB, 

(72) 

A^- (C, q) = [g^AA (.C, 9) + 9ab «\ 9)) - X+; 

(73) 

for free diblocks in a tube. 



Polymer-polymer and polymer-solvent correlations 

Substituting eqs (|i^ and ([S5t - ([55]) into expressions 
([70| and ([72]) we see that the functions A_ (0,0,(3^) are 
the same for free and tethered diblocks (bottlebrushes) in 
the case of a purely l ong itudinal ordering [n = 0, = 0) 
(see the well known [77| bold solid line 1 in Fig. 1 with a 
minimum at = = 3.785). It is worth to emphasize 
that the location of the minimum depends only on the 
length N of the side block (which scales the dimension- 
less quantity — q'^PN/6) and does depend neither on 
the grafting density nor on the value of the ^-parameter 
(temperature). This surprising similarity is due to the 
very nature of our approximation, which is valid only for 
strongly compressed bottlebrushes (k ^ 1). Indeed, in 
this case the entropic loss due to fixing one of the chain 
ends on the backbone is negligible as compared to that 
due to chain squeezing into a narrow tube, which is the 
same for free and tethered chains. For other concentra- 
tion modes having no angular but some radial modula- 
tion (77 = 0, s™ > 0) one has 

A_(0, Q') = A- (0, 0, Q2 + (s;r) VlG'^)) (74) 

(see definitions after eq l^^ . In other words, for these 
modes the instability curve A^{Q'^) — looks like that 
for the lamellar-like mode (tt, = 0, = 0) but shifted by 
(s;;")V(6k) (thin solid lines 1 and 2 in Fig. 1). 

Therefore, for the radially symmetric modes with 
{s^)^/{6Ql > K the function A_(0,C,Q^) reaches its 
minimum for Q = 0. For k, < kq = (s™)^/(6(5^ « 0.65, 
which is surely true within our approximation valid for 
K <C 1, the lamellar-like mode is the only one for which 
the function A_(0,s™,Q^) reaches its minimum at a fi- 
nite value of Q. Moreover, in this case the lamellar mode 
is dominant among those with no angular modulation. 
Indeed, as shown in Fig. 2, it is only the square averaged 
amplitude of this mode, which strongly increases when 



cl^Nx 



AB 



10.495. 



(75) 
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FIG. 1: The plots A_(Q^) for n = 0, to = 1 (the bold solid line), n = 0, to = 2 (thin solid lines) and n = 1, m = 1 
(the dashed lines). The labels 1, 2 and c correspond to the values oi k = 0.2, k ~ 0.05 and k = ki, respectfully. The 
bold dashed line corresponds to k = 0. The right figure shows the enhanced inset from the left one. The dotted 
straight lines are the guides for eyes to explain the meaning of the quantities Xc, Ql and Q^. 



The situation is rather different for the modes with 
angular modulation (n > 1). The relationship (fM| holds 
for these modes of free diblock copolymers in a tube also 
but for bottlebrush angular concentration modes one has 



A_(n,<\Q2)^A_(Q2 + (,™)2/^)^ 

A_(Q2) ^ (0Ar/^(QV2))"' - XAB. 



(76) 



In other words, the instability condition A_ = looks like 
that for the mixture of free blocks A and B but shifted 
by {s^Y /k (see the dashed lines in Fig. 1.). 

As described above, passing the instability threshold 
(|67l) results in formation of a lamellar-like or Janus- 
cylinders morphology depending on which function, (j70p 
or (1711) . has a deeper minimum. 

As shown in Fig. 3, if the shift (s™)^//c is small then 
the instability with respect to forming of Janus-cylinders 
appears first (with increase of xab)- With decrease 
of K the shift increases and we arrive at the value of 
K = Ki = (s})V(6Q?) « 0.0717 where both the lamellar 
and first angular modulated mode loose their stability si- 
multaneously at X = Xc With further decrease of n the 
lamellar phase becomes more thermodynamically stable. 
It follows from eq (fTTj) that n decreases when the distance 
a between the neighboring side chains of the bottlebrush 
increases (and, thus, grafting density decreases). Thus, 
the J anus- cylinders could be formed, in agreement with 
predictions of refs [1, IH, |4l[ in bottlebrushes with high 
grafting density. In bottlebrushes with lower grafting 
density we predict, contrary to refs d, H, 41 1, longitu- 
dinal lamellar-like (pearl-necklace) ordering rather than 



forming Janus-cylinder. 

It is worth to emphasize that our theory predicts a deep 
similarity between the polymer-polymer and polymer- 
solvent ordering in condensed bottlebrushes. Indeed, 
comparison of eqs ([701 and ([7T|) leads to conclusion that 
in symmetric bottlebrush the only difference in formation 
of the polymer concentration modulation described by 
the order parameter 4'„((r)) and solvent concentration 
modulation described by the order parameter ^_|_((r)) 
is energetic. More precisely, depending on whether the 
limit (iTSl) is reachable earlier or later than 



4>Nx^ 



Xc 



10.495, 



(77) 



the polymer-polymer or polymer-solvent inhomogeneities 
will occur but the type of morphology arising in both 
cases is expected to be the same. 



3. MOLECULAR DYNAMICS SIMULATIONS OF 
TWO-COMPONENT BOTTLE-BRUSHES 

Model and methods to analyze the results 

The strictly rigid backbone of the simulated bottle- 
brush polymer is simply taken as an immobile straight 
line in z-direction, where we also apply periodic bound- 
ary conditions, thus disregarding any end-effects. Using 
units of length such that the length parameter ctl j of the 
Lennard-Jones potential (see below) is taken to be unity, 
aLj = 1, we graft side chains (alternatingly of type A 
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FIG. 2: The square averaged amplitude of the lameUar mode (sohd) and first radiaUy (but not angular) modulated 
mode with n = 0, m = 2, s™ w 3.832 (dashed) for k = 1 (left), k — 0.5 (middle) and k = 0.25 (right). The first two 
cases are beyond validity of our approximation and are only suggestive. The values of A^{Q^) monotonously 
increase with increase of x which enables us not to label them. The bottom lines correspond to x = 0, the other 
lines do to X = Xc — 4'+^/10000, where i is the number of the curve from top to bottom. 




FIG. 3; The square averaged amplitude of the lamellar mode (solid) and first angular modulated mode with 
n — l,m — 1, 1.841 (dashed) for k — 0.11 (left), k = ki (middle) and k — 0.05 (right). For k — 0.1 the angular 

"Janus" mode appears (within the RPA) at x = Xi = fo {{sq)'^ /^2k) sa 8.471. The bottom lines correspond to 
X = 0, the other ones do to x = Xi - 4*+Vl0000 for k = 0.1 and x = Xc - 4*+7l0000 otherwise. 



and type B) at regular positions z^, = (ly — l)/cr, i' = 
1,2,..., 2M, with cr = 1.14 and a number 2M = 50 of 
grafted chains. The side chain length was chosen as N 
— 35. For comparison, also some results were obtained 
for a higher grafting density (cr — 1.51) and two lower 
ones (cr = 0.57 and a = 0.76); in the latter case, three 
chain lengths were studied {Na = Nb = N = 20, 35, and 
50, respectively). We are fully aware of the fact that our 
side chain lengths would not be long enough to study the 



asym ptot ic scaling behavior under good solvent condi- 



tions [47[, of course; however, for poor solvent conditions 
(as studied here) very large relaxation times would pre- 
vent us from reaching thermal equilibrium for longer side 
chains. Moreover, the range of N that we explore nicely 
corresponds to the range that can be probed in experi- 
ments 13, U, 



The side chains are modeled by the standard bead- 
spring model [23] already used in our previous work 
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5ll |55| . All beads interact with a truncated and shifted 
Lennard- Jones potential 



UlAt) = 



4eL,7[(cTLj/r)i2 - [aLj/rf] + C , r < 

, r > Vc 



(78) 

where Vc = 2.5aLj, and the constant C is defined such 
that ULjir = Tc) is continuous at the cutoff. Units are 
chosen such that ctlj = 1 irrespective of whether we deal 
with AA, BB or AB pairs, and also Boltzmann's constant 
fcs = 1, as well as the masses of the beads niA — "nis — 
m — 1. For the energy parameters, however, we choose 



AA 



1 



1/2; 3/4; 7/8; and 15/16, (79) 
1 the problem 



respectively. Remember that for e;^^ 
reduces to a single-component brush. 

As usual [24, 51, 55] the connectivity of the beads along 
a side chain is maintained by the "finitely extensible non- 
linear elastic" (FENE) potential, 

UFENE{r) = -hrl ln[l - [r/ro f], < r < ro , (80) 

with the standard choice of parameters, ro = 1.5, and 
fc = 30, respectively. Since UpENEir > ro) = oo the 
bond lengths are constrained. 

Simulations were carried out controlling temperature 
by the Langevin thermostat 24 , 5ll [5^ . The equation of 
motion of a bead then reads 



mj—+r.{t), 



(81) 



where t denotes the time, Ui is the total potential acting 
on the i-th bead (which is at the position r^), 7 is the 
friction coefficient, and Ti{t) the random force. As is 
well known, 7 and Ti{t) are related by the fluctuation- 
dissipation relation 



(f,(t) • fj{t')) = 6kBT-fd^jd{t' - t). 



(82) 



Eqs. I81I82I are numerically integrated using the GRO- 
MACS package, applying the leap frog algorithm 87. SSf. 
Following previous work 2J, |5l|, |55|, 7 = 0.5 was cho- 
sen, and a time step At = 0.006t where the Molec- 
ular Dynamics(MD) time unit was chosen as r = 
«,/eij)V2 = l. 

As was briefly mentioned above, equilibration of our 
model system is rather difficult. Typically we started by 
equilibrating the system first at a temperature T — 3.0, 
which is known [89| to correspond roughly to the Theta 
temperature of homopolymers in dilute solution. 500 dif- 
ferent configurations at T = 3.0 were equilibrated using 
runs extending over a time range of 30 x lO^r. All these 



configurations are then cooled down from T = 3.0 to 
T = 1.5 with temperature steps AT — 0.1. At each 
temperature 2 x 10^ steps were used for further equi- 
libration. As an example. Figs. 4 and 5 show typical 
snapshot pictures of the simulated bottle brushes. At 
T = 3.0 the individual side chains take irregular confor- 
mations, stretching out in the directions perpendicular to 
the backbone. Obviously, these conformations are rather 
disordered, the contacts between monomers of different 
chains are too scarce to induce significant microphase 
separation, and even on local scales the composition of 
the bottle brush is rather randomly mixed. At T = 1.5 
the situation is different, however: for a = 0.57 one can 
clearly see the presence of a pearl-necklace type structure, 
as far as the total density is concerned; typically these 
pearls sometimes are almost exclusively formed from one 
type of chain only, implying a periodic concentration vari- 
ation along the z-axis. For a = 0.76, however, the inter- 
pretation of many snapshot pictures rather suggest that 
"Janus pearls" (containing an A- rich and a B-rich part, 
separated by an interface) have formed. For a = 1.14, 
already rather elongated A-rich and B-rich parts (remi- 
niscent of "Janus cylinders") are visible, but there still 
occur strong undulations in the local monomer density in 
z-direction. These undulations for a = 1.51 clearly are 
somewhat weaker, although the density variation still is 
far from a uniformly filled cylinder. While some snap- 
shots suggest that a "Janus cylinder" type ordering per- 
sists over the full lengths of the studied backbone length, 
we emphasize the caveat that in Figs. 4, 5 one cannot 
really distinguish whether a local cross section perpen- 
dicular to the cylinder axis resembles more a sphere (as 
appropriate for a cylinder) or a dumbbell, for instance. 

In order to provide a quantitative analysis of the vari- 
ous types of short range order we define correlation func- 
tions adapted to the cylindrical symmetry of our prob- 
lem, by suitable generalization of the two-point correla- 
tion function in fluid systems. 



— Wo ^fl 

GMSzJ^) = -^(E E S(z. - - Sz) (83) 

Sicj),-cbj-S(b)), a,l5 = A,B 

In Eq. [M] the first sum runs over all N ^ particles of 
type a, which have the cylindrical coordinates Zi,ri,4>i, 
while the second some runs over all particles of 
type which have the cylindrical coordinates Zj, rj, 
Na = = NM. Fig. 6 shows this correlation func- 
tion GaaISz, S(/)) for a typical case: as expected, in 
the temperature region close to the 0-temperature this 
correlation decays to its asymptotic value for large Sz, 
Gaf3{Sz — >■ 00, 6(j)) = 1 if Sz > 10, while a nontrivial an- 
gular dependence occurs for smaller Sz. Basically, this 
structure is due to intra-chain-correlations. The corre- 
lation function GAsiSz, S(j)) in the region Sz < 10 is 
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predominantly smaller than unity, the chains avoid un- 
favorable contacts already at this rather high temper- 
ature. However, the angular dependence occurring at 
T = 2.9 in GAB{^z,5(j)) for 5z > 10 clearly is due to a 
lack of statistical accuracy, as a comparison of data for 
individual runs shows. Similarly, the weaker deviations 
of Gaa{Sz > 15, S(j)) from unity should be disregarded 
for the same reason either. While at low temperatures 
(T = 1.5, see Fig. 6 clearly a more systematic correla- 
tion over a larger range of Sz develops, the systematic 
finite size effects due to the periodic boundary condition 
in 0-direction become a concern, since 

{Sz, S^) = G^p (L, - 5z, 5(t>) (84) 

where the length Lf, of the backbone is Lb — 2M /a (in 
the examples shown in Fig. 6 we have Lb = 87.72, i.e., 
Gap{Sz, 54)) is shown for < < -Lb/2). 

The conclusion of the discussion of Fig. 6 is that a 
substantially larger statistical effort would be needed to 
sample Gaf3{Sz, Scj)) with meaningful accuracy. Therefore 
it was decided to focus on a more meaningful average 
information, defined as 



on gap{Sz,m), we also consider the structure factor 
Sap{q) of the bottle brushes; as is well-known, for 
one-component bottle-brushes ample information on the 
structure factor S{q) is available from scattering exper- 
iments [13, [lH, and for two-component bottle-brushes 
relevant information could be obtained from scattering 
experiments as well if one species (A or B) is a deuterated 
polymer. Assuming that the wave vector of the scatter- 
ing q is also oriented along the z-axis of the bottle-brush 
polymer, the partial structure factor Sa^iq) becomes 

k=l 1=1 

where a/3 G A,B, N — Na + Nb- Again it is useful 
to transform to structure factors that relate to total den- 
sity fluctuations {Snn{q)) and concentration fluctuations 
(Scciq)), respectively, 

Snniq) = SAAiq) + ^SAB^q) + SBsiq), (89) 



Gafs{Sz,ni) = J Gai3{Sz,S(j)) cos{mS(f>)d(f>,m — 0.1,2 



(85) 

The choice m = simply means that the angular dif- 
ference dcj) between the coordinates r'i,rj is not at all 
taken into consideration. The choice m = 1 is used in 
the expectation that this correlation then will yield some 
information whether or not Janus-cylinder type ordering 
occurs, while the choice m — 2 could be useful, if "Janus- 
dumbbell" and double-cylinder-like structures occur (see 
also the discussion in Sec. 2.2). 

We have also used these correlations to transform to 
correlation functions that reflect fluctuations of total 
number density (gnn{Sz,m)) and relative concentration 
{gcc{Sz, m)). Following Bhatia and Thornton [q^I we ob- 
tain 

gnn{Sz,m)^ ^ XaXf)Gal3{Sz,m) (86) 
a,l3=A,B 

and 

gcc{Sz,m)^ ^ yayi3Gafi{5z,m) (87) 

a,p=A.,B 

where in our case xa = xb — 1/2 and yA — ^yB — 1/2. 
Of course, gnn(5z, m) can be computed directly, consider- 
ing a correlation Gai3 where the sums over i and j fEa. l84p 
both run over all particles, irrespective of whether they 
are of type A or B. 

Since it is not clear whether measurements will 
become possible that will yield information directly 



Scc{q) = x\SAA{q) + x\SBB{q) - 2xAXBSAB{q)- (90) 

In an actual experiment, of course, it would be more 
natural to consider the situation when the orientation 
of the scattering vector q is fixed by the experimental 
setup, while the orientations of the rigid backbones of the 
bottlcbrush polymers in solutions are randomly oriented. 
Then one would need to consider a structure factor 

Sc.p{q) = = ^ 5](cxp[ig • (ffe - n)]) , (91) 

k=l 1=1 

where it is understood that the average (...) includes 
an average over the orientation of q. This struc- 
ture factor Sap{q) can be interpreted as the Fourier 
transform of a correlation function gQ,^(Ar), where 
Ar is the absolute value of the distance between two 
sites of monomers ri,rj in the bottle brush, Ar = 
{xj — XiY + {yj — yiY + (^i ~ ^i)^- This correlation 
(7a^(Ar) was presented and discussed already in our pre- 
liminary communication [ssf . and hence will not be dis- 
cussed further here. We only note that for large val- 
ues of Ar we essentially must have Ar « \zj — Zi\ and 
then 5Q/3(Ar) « gap{5z = Ar), since the differences 
Xj — Xi , yj — yi must be relatively much smaller (of the 
order of the cross-sectional radius R of the cylindrical 
brush). However, for small Ar the behavior of cjapi^r) 
is needed to estimate the number of contacts n^^^ which 
is defined as 
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An 



(92) 



Eq. [92] means that a pair of particles (a,/3) is defined to 
have a pairwise contact if their distance is less than rcon- 
We have followed the criterion due to Stillinger to 

take Tcon = ^-^ffLJ- 



RESULTS 

We first focus on the dependence of the number of 
contacts = rt^^ and on temperature T and 

interaction energy cab (remember e^A — ^bb = e = 1)i 
Fig. 7, to clarify what is the range of interest of those pa- 
rameters, where the microphase separation occurs. We 
recognize that for T > 3 all contact numbers start to be- 
come independent of T when T is increased. This behav- 
ior is easily interpreted by the fact that for T > Q(B « 3) 
the repulsive interactions among the beads dominate, the 
side chains essentially take on self-avoiding walk-like con- 
figurations. Of course, each interior monomer of a side 
chain has at least two monomers of the same type (A 
or B) which are nearest neighbors along the same side 
chains, so we expect that = cannot be less than 
two. Fig. 7 shows that the actual numbers — 
at temperatures slightly above T = are about twice as 
large, n"^ « 4. In contrast, the numbers of heterocon- 
tacts between chains of different type are at least an order 
of magnitude smaller, from n^^(T = 3.3, €a = \) ~ 0.20 
to n^^JT = 3.3, eAB = 15/16) w 0.42. It is interesting to 
note that the number = nfj^ increases strongly with 
decreasing temperature for T < 2.0, which we identify as 
the regime of temperatures where for the chosen (short!) 
chain length of the side chains a rather dense and ap- 
proximately cylindrical structure of the bottle-brush can 
be identified. Note that = n^J^ at these low tem- 
peratures decreases systematically with increasing cab, 
while increases with €ab- However, the dependence 



A A 

con 



,AB 



) on eAB is 



of the total number of contacts (n^ 
much weaker. 

On the other hand, the temperature dependence of 
'^con depends very strongly on cab ■ while for 
CAB < 3/4 stays almost independent of T, even for low T 
where the brush is in a rather dense state, for cab = 7/8 
and EAB — 15/16 we find a strong increase of nfj^ with 
decreasing T. Tentatively, this increase can be attributed 
to the formation of AB-interfaces in the collapsed parts 
of the bottle-brush. 

We now turn to the behavior of the correlation func- 
tions g„n((5z, m) and gcci^z, m) and focus on the grafting 
density a — 1.14 and compare two choices of cab, ^ab = 
3/4 and tAB = 15/16 (Figs. 8-11). We see that for 
^AB = 3/4 (Fig. 8) in the average density fluctuation 



at high T there is a rapid decay of gnn{Sz,m = 0) to 
unity (deviations from unity of order 0.01 are just statis- 
tical noise) , while at lower temperature a clear minimum 
near 5z = 10 followed by a maximum near 5z = 25 — 30 
develops, indicating rather long wavelength fluctuations 
in the thickness of the collapsed cylinder. This long- 
wavelength structure is not present in gnn{Sz,m — 1), 
Fig. 8b, while gnn{5z,m — 2), Fig. 8c, indicates that a 
uniform correlation gradually develops, which does not 
occur for cab = 15/16, however (Fig. 10). While in the 
average concentration correlation gcc{5z, m = 0) one sees 
a long wavelength periodicity for cab =3/4 (Fig. 9a), 
as well as for tAB = 15/16, (Fig. 11a), this periodicity 
is basically independent of temperature and its ampli- 
tude is very small (of order 0.005). Thus, this variation 
clearly must be discarded as being due to somewhat in- 
sufficient preparation of the initial states. In contrast. 
Fig. 9b shows that gcci^z, 1) develops much stronger or- 
der of uniform sign, indicating that this system shows 
a tendency to develop Janus-cylinder type short range 
order. While at €ab — 7/8 (not shown) this tendency 
is also present, though less developed, for €ab — 15/16 
(fig. lib) the ordering tendency is much weaker: this 
indicates that the effective x-parameter for phase sepa- 
ration between A and B at eab = 15/16 and = 35 
may be too weak to cause micro-phase separation. The 
correlation function gcc{Sz,m — 2), Figs. 9c and 11c 
show a rapid decay to zero in both cases, independent 
of temperature, indicative of the fact that this correla- 
tion measures only some intra-chain correlation effects 
in these cases, and no collective effects are detected by 
gcc{Sz, m = 2) at all. 

When we study still larger grafting densities, such as 
a = 1.51 (Figs. 12, 13), these conclusions are confirmed: 
there is no significant long wave length correlation in ei- 
ther uniform density or concentration, but a pronounced 
tendency to form Janus-cylinder like short-range order 
develops, together with a long range correlation in the 
density, seen in (7„„((5z,m = 2), which is due to a cross- 
section of the cylindrical brush which is locally non- 
spherical. The obvious interpretation of Fig. 8c) and 
12c) is that the local Janus- type order is not strictly of 
Janus cylinder type but rather of Janus dumbbell type, 
for EAB = 3/4, as it was suggested already in our pre- 
liminary communication 55| where other quantities for 
the case sab = 1/2 were studied. This conclusion is 
corroborated by the concentration correlations (Fig. 13), 
which reveal neither a significant long range correlation in 
gcci^z, 171 = 0) nor in gcdSz, m > 2) while gcd^z, m = 1) 
exhibits the development of very pronounced Janus-type 
order. 

For considerably lower grafting density, such as cr = 
0.57 (not shown) or cr = 0.76 (Figs. 14,15) the behav- 
ior is very different, however: both f;„„((5z,rn = 0) and 
gcc{Sz, m = 0) do not exhibit much structure at high 
temperature, but develop an oscillation (minimum near 
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5z = 10, maximum near 6z = 20) when T is lowered to- 
wards T = 1.5. For longer chain lengths (N = 50, not 
shown) the behavior is qualitatively similar, but the min- 
imum now occurs near Sz = 12 and the maximum near 
5z = 24, and for N = 20 (not shown) gcc{5z,m — 0) 
has its minimum at 5z = 6, the maximum at 6z — 12. 
Thus, there is a clear dependence on this weak correla- 
tion seen in Figs. 14a, 15a on the length of the side- 
chains. In g.mi{Sz,m = 1) there occurs only a fast decay 
to zero for N=20 (not shown), while for N=35 a mini- 
mum develops for 5z = 10 (Fig. 14b) and for N=50 this 
minimum occurs for dz = 12, i.e., these structures are 
strongly correlated with the structure that develops for 
gnniSz,m = 0). Taken together with the observation 
(Fig. 5b) that the chains form a pearl-necklace struc- 
ture, we conclude that typically the "mass" in the pearls 
is not distributed spherically symmetric around the back- 
bone (the z-axis) . This feature is also suggested by visual 
inspection of the snapshots (Fig. 5a, 5b), but Fig. 14 im- 
plies that it is not an accidental observation of some par- 
ticular snapshots, but a statistically significant physical 
effect. 

Turning to the concentration correlation gcc{Sz,m = 
1) we notice that for high temperatures a rapid decay to 
zero occurs but this decay becomes significantly slower as 
T is lowered (Fig. 15b). Again the comparison with the 
corresponding data for TV = 20 (not shown) and N=50 
(not shown) reveals a qualitatively similar picture (for N 
= 50 the decay is so slow, however, that our choice of 2M 
= 50 would not be large enough to avoid artefacts due 
to the periodic boundary conditions) . Turning finally at- 
tention to gnn{5z,m = 2) and gcc{,5z,m — 2), Figs. 14c 
and 15c, we notice a monotonous decay to zero, which be- 
comes somewhat slower when T is decreased. For N = 20 
(not shown) these correlations are almost independent of 
T, decay to zero for 5z > 5, while for N = 50 (not shown) 
gnni^z, m = 2) develops a slow decay (which is absent in 
gcc{Sz,m — 2), however). Clearly these data are deli- 
cate to interpret, since upon lowering the temperature 
for this larger chain length one presumably crosses the 
boundary from an essentially uniform cylindrical struc- 
ture (as function o f z) at higher temperatures to rather 
elongated clusters [56|. If one compares the behavior to 
the case a = 0.57 (here some snapshot pictures. Fig. 
5a, suggest a clear alternation of A-rich and B-rich re- 
gions along the z-axis), one does not find any significant 
structure of longer range in either gnniSz,m = 1,2) or 
gcc{Sz, m — 1,2) [in order to save space, the correspond- 
ing data are not shown], and only both gmi{dz,m = 0) 
and gcciSz,m = 0) exhibit the development of periodic 
structures, similar to what is seen in Figs. 14a, 15a. 

Finally Fig. 16 shows a few typical results of structure 
factors Scc{q)- Note that due to the periodic bound- 
ary conditions q is "quantized" in multiples of 27r/L6 = 
TTa/M, with 2M = 50 in Fig. 16 (a,b) while 2M = 100 
in Fig. 16 (c,d). 



We find that the position of the main peak occurs at 
Qmax ~ 0.3, (which corresponds to Qmax * Rq = 1.92, 
where Rq would be the gyration radius of a diblock 
macromolecule Aj^Bn with N = 35, as used here, at 
the Theta point, T — 3.0) in all cases, it is only the 
peak intensity that strongly grows with increasing a or 
changing the temperature. Such a finding is quite non- 
trivial, and without guidance by the theoretical treatment 
of the previous section,would have been somewhat unex- 
pected. Moreover, it is worth to notice that the theoreti- 
cal and simulation values for x = (qmax * Ro)^ (3.785 and 
1.922 w 3.685, respectively) are in remarkable agreement. 
Since the theory is done for temperatures slightly below 
the Theta temperature but very long side chain lengths, 
while the simulation is done for temperatures down to one 
half the Theta temperature but rather short side chains, 
the agreement is perhaps better than one would have ex- 
pected. 

In previous work on cylindrical bottle-brushes with a 
sing le type of side chain under poor solvent conditions 
[56| it was found that for 0.4 < cr < 0.8 clusters elon- 
gated in z-direction occur, which have a gyration radius 
of order (i?^^)^^^ ~ 10. These elongated clusters are also 
visible in the snapshot pictures (Figs. 5a, b). For larger 
a the range over which the concentration distribution is 
uniform (as measured by gcc{5z,m = 1) see Fig. 9b) in- 
volves a similar length scale. However, the lack of a more 
pronounced temperature dependence of Sccil) is disturb- 
ing, and lacks an explanation. We conclude that in any 
case it seems difficult to clarify the nature of short range 
order in binary bottle-brush polymers from an analysis 
of Scc{q)- 

In order to get a more precise picture of the mesophase 
ordering in reciprocal space from the simulations, which 
could be compared to the corresponding results of the 
analytical calculations of Sec. 2 (e.g.. Figs. 2, 3), it is 
useful to carry out numerical Fourier transforms of the 
data on the correlations in real space, presented in Figs. 8 
-15. For example. Fig. 17 shows the Fourier transforms 
Snnilj 1) ^-nd Snniq, 2) which are the Fourier transforms 
of gnni^z,!) and f/„„((5z,2), for a representative selec- 
tion of our data, while Fig. 18 shows the Fourier trans- 
form of the corresponding concentration fluctuation cor- 
relations {ScciQ, 1) is the Fourier transform of gcciSz, 1), 
while Scc{q, 2) is the Fourier transform of gcc{Sz, 2)}. One 
can see that for a — 0.76 the density fluctuation struc- 
ture factors Snniq,!) and ^^^(g, 2) show a very pro- 
nounced increase at small q as the temperature is low- 
ered, while for a = 1.14 this increase is much weaker, 
and for Sn7i{q,l) it is evident that a small peak grows 
at nonzero q only (g « 0.3), compatible with the qmax 
extracted from Snn{q) as shown in Fig. 16. We conclude 
from these data that density variations play a central 
role for mesophase formation for a = 0.76, but less so for 
a = 1.14. The concentration fluctuation structure factor 
5cc('Z, 0) exhibits for a = 0.76 a strongly growing peak 
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near qm ~ 0.3 (not shown since the data are quahtatively 
similar to Fig. 16b), while Scdq,^) for a = 1.14 exhibits 
a very weak temperature dependence only (not shown), 
as it must be, considering that gcc{5z^Q) in this case is 
basically independent of temperature (Fig. 9). However, 
the trend seen for Scc{q^ 1), which exhibits a growth for 
small q with decreasing temperature for all the cases in- 
cluded in Fig. 18, indicates a tendency for local Janus 
cylinder type order in these cases, while long range order 
clearly is not established. We have not included Scc{q, 2) 
into this study, since the data for gcc{Sz, 2) {Fig. 9c, 11c, 
13c} indicated already hardly any interesting tempera- 
ture dependence to be present. 

The data shown in Figs. 16 - 18 indicate that in the case 
a = 0.76 density takes the leading role in the mesophase 
ordering of the cylindrical brush, and the concentration 
ordering adjusts to it, while for a = 1.14 concentration 
variations along the backbone drive the mesophase or- 
dering, with a strong coupling between density variations 
and concentration variations still being present. In this 
case, a periodic modulation of density along the z-axis 
still grows when the temperature is lowered (Fig. 17a), 
and correspondingly also Scdq) has a pronounced peak 
near q q^ ^ 0.3 (Fig. 16 c,d). On the other hand, for 
a = 1.51 a rather different behavior is found (Fig. 19): 
while for SnniqA) only a weak peak at about q sa 0.2 
grows when T is lowered (Fig. 19a), the growth of both 
S'„„(g, 2) and ScdqA) for q — > is very strong. We 
interpret these observations as a tendency towards the 
formation of Janus dumbbell- type structures. 



4. CONCLUSIONS 

In this paper the microphase separation in binary 
(A,B) bottlebrush polymers with rigid backbones driven 
by decrease of temperature (i.e., variation of the sol- 
vent quality from Theta solvents to poor solvents) was 
studied, considering also the variation of grafting density 
along the backbone of these cylindrical brushes. Two 
complementary theoretical methods were used, namely, 

(i) a suitable extension of the random phase approxi- 
mation (RPA), and 

(ii) Molecular Dynamics (MD) simulations. 

The first approach has the merit that arbitrarily long 
chain lengths of the side chains can be considered, as well 
as infinitely long backbones (so that end effects or effects 
of boundary conditions along the backbone do not mat- 
ter). Besides, the approach results in some explicit ex- 
pressions for the correlation functions under study. How- 
ever, the disadvantage of the approach is that it is basi- 
cally a linear stability analysis around the homogeneous 
state, strongly nonlinear effects are out of consideration. 



The MD approach, on the other hand, in principle 
can take both nontrivial correlations and nonlinear ef- 
fects into account, but in practice (to avoid excessive re- 
quests for computer time resources) is limited to rather 
short length of the side chains, and also the length of the 
(rigid) backbone is finite (end effects then are avoided 
by a periodic boundary conditions, but the price to be 
paid is that reciprocal space in the z-direction along the 
backbone is discretized) . Also, non- negligible statistical 
errors (and systematic errors due to insufficient length 
of the MD trajectories, which affect the establishment 
of full thermal equilibrium) restrict the accuracy of the 
results that can be gotten. 

Whereas the RPA treatment assumes that the 
monomer density in the radial direction perpendicular 
to the cylinder axis is uniform and constant (up to the 
cylinder radius, while it is zero outside of the cylinder 
surface), in the model studied by MD this assumption 
clearly does not apply, as has been shown in our previous 
work [51I [13, [HBl- Rather, it was found that there occur 
density oscillations very close to the cylinder axis (rem- 
iniscent of the layering of fluid particles adjacent to the 
hard wall) , and for the chain lengths one can study only 
for a rather small regime of radial distances the density 
profile then is approximately constant. One approaches 
rather fast the regime where the density continuously 
drops to zero, similar to the density profiles across inter- 
faces in phase-separated polymer solvent systems, where 
the density decreases from its value in the melt towards 
zero over a distance of several monomer diameters. This 
low temperature behavior in the simulation comes clos- 
est to the situation assumed in the RPA calculations. At 
higher temperatures, close to the O temperature, the in- 
terfacial profile is so broadened that a flat region where 
the density has reached the melt density no longer can 
be identified. Given the different temperature and chain 
length regimes for which theory and simulation access 
the microphase separation conditions, it is very gratify- 
ing that there is not only a good qualitative agreement 
between both, but furthermore the wavelength of insta- 
bility towards the microphase separation agrees quanti- 
tatively between both. 

Qualitatively, the analytical theory describes similar 
trends as seen in the simulation: there is a competition 
between microphase separation developing periodic order 
in axial direction, i.e., the system tends to develop long 
range order with a particular periodicity (the wavelength 
of this periodic modulation being controlled by the gyra- 
tion radius of the side chains, while the strength of the 
modulation depends on temperature and grafting den- 
sity, as the analytical theory compellingly shows, and as 
is confirmed by our simulations) and microphase separa- 
tion of the Janus-type, characterized by g = (the cross 
section of the dumbbell, depending on the ratio eAs/^-, 
in the strongly segregated regime). 
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The theory shows that the instabihty leading to the 
latter type of order prevails if the grafting density is suf- 
ficiently high (and/or the side chain length sufficiently 
long). Both the snapshot pictures of the simulated 
bottle-brushes (Figs. 4,5) and the quantitative analysis 
of the simulations in terms of the correlation functions 
of density fluctuations g„„((5z,m) [with m = 0,1,2] and 
concentration fluctuations (or their Fourier transforms) 
are compatible with such an interpretation. 

We emphasize that both the simulation and the ana- 
lytical theory have considered a particularly symmetric 
case, where both types of side chains have equal length 
{Na = Nb) and also their number was taken exactly the 
same (and moreover the same solvent quality was chosen 
for both types of chains, caa = ^bb)- Clearly, when 
one would like to discuss particular experimental sys- 
tems, these rather special conditions need to be relaxed. 
Also, we have restricted attention to the case where the 
side chains arc perfectly flexible while the backbone was 
assumed to be rigid. 

For many real bottle-brush systems it is probably 
more realistic to assume backbones that are semi-flexible 
rather than completely rigid, and also the local intrinsic 
stiffness of the side chains may play a role. Thus, the 
treatment presented in our paper can be taken as a flrst 
step towards a more complete description of less complex 
but interesting systems only. 
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FIG. 4: (Color online) Snapshot picture of bottle-brush 
polymers at cr = 1.14, N = 35, eAB = 3/4 for T = 3.0 

(upper part) and T = 1.5 (lower part). A and B 
monomers are distinguished by different color (or light 
grey vs. dark grey, respectively) 
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(b) 
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(d) 

FIG. 5: (Color online) Snapshot picture of bottle-brush 
polymers at eAB — 1/2, N = 35, T = 1.5 and four 
different values of ct : cr = 0.57 (a), 0.76 (b), 1.14 (c) 
and 1.51 (d) 
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FIG. 7: (Color online) Number of contacts n 
(a) and (b) plotted vs temperature, for a 
and several choices of £^15,6^5 = 1/2, 3/4, 7/8 and 
15/16, as indicated. For comparison, the case a = 1.51 
^AB =3/4 also is included. All data refer to TV = 35 

and 2M = 100 was chosen (in order to minimize 
artefacts due to the periodic boundary conditions). 
Note that in part (b) the numbers of contacts are 
multiplied by a factor 10 throughout. 
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FIG. 6: (Color online) Correlation function 
Gaa{^z, S(j)) plotted vs. Sz for the case a — 0.57, N = 
50, eab — 1/2, and the temperatures T — 2.9 (a) and 
1.5 (b). Several choices of the angle 5(f> — kir/SG, 
/c = 1, 4, 8, 12, 16, 20, 24, 28, 32 and 36 are shown, as 
indicated (c). Correlation function Gab{5z,5<1)) plotted 
vs. 5z for the case a — 0.57, N = 50, eab = 1/2 and 

r = 2.9 . 
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FIG. 8: (Color online) Correlation functions (7„„((5z, m) 
relating to density fluctuations plotted vs. 6z and the 
case €AB — 3/4, N = 35, cr = 1.14 for m = (a), m — 1 

(b) and m — 2 (c). Eight temperatures 
T = 1.5, 1.8, 2.1, 2.4, 2.7, 3.0, 3.3 and 4.0 are included, as 
indicated. 
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FIG. 9: (Color online) Same as Fig. 8, but for the 
correlation function gcciSz,m) relating to concentration 
fluctuations. 
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FIG. 10: (Color online) Same as Fig. 8, but for the 
case Cab = 15/16. 
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FIG. 11: (Color online) Same as Fig. 9, but for the 
case eAB = 15/16. 
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FIG. 12: (Color online) Same as Fig. 8, but for the 
case €ab = 3/4, N — 35,a = 1.51. Six temperatures 
from T — 1.5 to T = 3.0 are shown, as indicated. 
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FIG. 13: (Color online) Same as Fig. 9, but for the 
case eab — 3/4, = 35, a = 1.51. 
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FIG. 14: (Color online) Same as Fig. 8, but for the 
case eAB = 1/2, iV = 35, = 0.76. 
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FIG. 15: (Color online) Same as Fig. 9, but for the 
case eAB = 1/2, iV = 35, = 0.76. 
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FIG. 16: (Color online) Structure factor Scc{q) plotted 
vs. qioT a = 0.57, N = 35, e^s = 1/2 (a), a = 0.76, N 
= 35, CAB = 1/2 (b) a ^ 1.14, N = 35, e^B = 1/2 (c), 
and a = 1.14, N = 35, e^s 7/8 (d). The positions of 
the main peaks occur near qmax ~ 0.3 in all cases, 
corresponding to a characteristic wavelength 

A « 27r/qmax ~ 21. 
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FIG. 17: (Color online) Fourier transforms of density 

fluctuation correlation functions, Snnil, 1) and 
Snn{q, 2), as defined in the text, for the case N = 35, 
a = 1.14, EAB = 3/4 (a, b) and for the case N = 35, 
a ^ 0.76, eAB = 1/2 (c, d). 
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FIG. 18: (Color online) Fourier transforms of 
concentration fluctuation correlation functions Scdq, 1) 
for the cases iV = 35, cr 1.14, e^s = 3/4 (a), N = 35, 
a = 1.14, EAB = 15/16 (b) and TV = 35, ct = 0.76, 
eAB = 1/2 (c). 
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FIG. 19: (Color online) Fourier transforms of density 

fluctuation correlation functions S'„„(g, 1) (a) and 
S'„„(g, 2) (b), and concentration fluctuation correlation 
functions Scc{q, 1) (c) for the case N — 35, a — 1.51, 
,and eAB = 3/4. 



25 



M. Zhang and A. H. E. Miiller, Polym. Sci. Part. A: 
Polyni. Chcm. 43, 3461 (2005) 

S. S. Sheiko, B. S. Sumerlin, and K. Matyjaszewski, 
Progr. Polym. Sci. 33, 759 (2008) 

M. Wintermantel, M. Schmidt, Y. Tsukahara, K. Kaji- 
wara and S. Kahijiya, Makromol. Chem., Rapid Com- 
mun. 15, 279 (1994) 

M. Wintermantel, M. Gerlc, K. Fischer, M. Schmidt, I. 
Wataoka, H. Urakawa, K. Kajiware, and Y. Tsukahara, 
Macromolecules 29, 978 (1996) 

K. L. Beers, S. G. Gaynor, K. Matyaszewski, S. S. Sheiko, 

and M. MoUer, Macromolecules 31, 9431 (1998) 
S. Lccommandoux, F. Chilecot, R. Borsali, M. Schap- 
paclier, A. Deflieux, A. Brulet and J. P. Cotton, Macro- 
molecules 35, 8878 (2002) 

T. Stephan, S. Muth, and M. Schmidt, Macromolecules 

35, 9875 (2002) 

Y. Liu, V. Abetz and A. H. E. Miiller, Macromolecules 

36, 7894 (2003) 

C. Li, N. Gunari, K. Fischer, A. Janshoff, and M. 
Schmidt, Angew. Chem. Int. Ed. 43, 1101 (2004) 
S. Rathgeber, T. Pakula, K. Matyaszewski, and K. L. 
Beers, J. Chem. Phys. 122, 129404 (2005) 

B. Zhang, F. Grohn, ,J. S. Pedersen, K. Fischer and M. 
Schmidt, Macromolecules 39, 8440 (2006) 

L. Feuz, P. Strunz, T. Gene, M. Textor and O. Borisov, 
Eur. Phys. J. E23, 237 (2007) 

H. I. Lee, K. Matyaszewski, S. Yu-Su, and S. S. Sheiko, 
Macromolecules 41, 6073 (2008) 

R. V. lozzo (ed.) Proteoglycans: Structure, Biology and 

Molecular Interactions (Marcel Dekker, New York, 2000) 

H. Muir, Biocheni. Soc. Trans. 11, 613 (1983) 

N. C. Kaneider, S. Dunzendorfer, and C. J. Wiedermann, 

Biochemistry 43, 237 (2004) 

J. Klein, Science 323, 47 (2009) 

T. M. Birshtein and E. B. ZhuUna, Polymer 25, 1453 

(1984) 

T. Witten and P. A. Pincus, Macromolecules 19, 2509 

(1986) 

T. M. Birshtein, O. V. Borisov, E. B. Zhulina A. R. 
Khokhlov, and T. A. Yurasowa, Polym. Sci. USSR 29, 
1293 (1987) 

Z. G. Wang and S. A. Safran, J. Chem. Phys. 89, 5323 

(1988) 

C. Ligoure and L. Leibler, Macromoelcules 23, 5044 
(1990) 

R. C. Ball, J. F. Marko, S. T. Milner, and T. A. Witten, 
Macromolecules 24, 693 (1991) 

M. Murat and G. S. Grest, Macromolecules 24, 704 
(1991) 

N. Dan and M. Tirrell, Macromolecules 25, 2980 (1992) 
C. M. Wijmans and E. B. Zhulina, Macromolecules 26, 
7214 (1993) 

G. H. Fredricksen, Macromolecules 26, 2825 (1993) 

H. Li and T. A. Witten, Macromolecules 27, 449 (1994) 
E. B. Zhuhna and T. A. Vilgis, Macromolecules 29, 2605 
(1996) 

Y. Rouault and O. V. Borisov, Macromolecules 29, 2605 
(1996) 

E. M. Sevick, Macromolecules 29, 6952 (1996) 

M. Saariaho, O. Ikkala, I. Szleifer, I. Erukhimovich, and 



G. ten Brinke, ,J. Chem. Phys. 107, 3267 (1997) 
[33] M. Saariaho, I. Szleifer, O. Ikkala and G. ten Brinke, 

Macromol. Theory Simul. 7, 211 (1998) 
[34] Y. Rouault, Macromol. Theory Simul, 7, 359 (1998) 
[35] M. Saariaho, A. Subbotin, I. Szleifer, O. Ikkala, and G. 

ten Brinke, Macromolecules 32, 4439 (1999) 
[36] K. Shiokawa, K. Itoh and N. Nemoto, J. Chem. Phys. 

Ill, 8165 (1999) 
[37] A. Subbotin, M. Saariaho, O. Ikkala and G. ten Brinke, 

Macromolecules 33, 3447, (2000) 
[38] P. G. Khalatur, D. G. Shirvanyanz, N. Y. Staroviotovo, 

and A. R. Khokhlov, Macromol. Theory Simul. 9, 141 

(2000) 

[39] R. Stepanyan, A. Subbotin, and G. ten Brinke, Macro- 
molecules 35, 5640 (2002) 
[40] N. A. Denesyuk, Phys. Rev. E67, 051803 (2003) 
[41] J. de Jong and G. ten Brinke, Macromol. Theory Simul. 
13, 318 (2004) 

[42] S. EUi, F. GanazzoU, E. G. Timoshenko, Y. A. 
Kuznetsov, and R. Connolly, J. Chem. Phys. 120, 6257 

(2004) 

[43] S. S. Sheiko, O. V. Borisov, S. A. Prokhorova, and M. 

Moller, Eur. Phys. J. E 13, 125 (2004) 
[44] R. Connolly, G. Bellesia, E. G. Timoshenko, Y. A. 

Kuznetsov, S. Elli, and G. Ganazzoli, Macromolecules 

38, 5288 (2005) 
[45] A. Yethiraj, J. Chem. Phys. 125, 204901 (2006) 
[46] H.-P. Hsu, W. Paul, and K. Binder, Europhys. Lett. 76, 

526 (2006) 

[47] H.-P. Hsu, W. Paul, and K. Binder, Macromol. Theory 

Simul. 16, 660 (2007) 
[48] A. V. Subbotin and A. N. Semenov, Polym. Sci., Ser. 

A49, 1326 (2007) 
[49] H.-P. Hsu, W. Paul, and K. Binder, J. Chem. Phys. 129, 

204904 (2007) 

[50] I. I. Potemkin and V. V. Palyulin, Polym. Sci. Ser. A51, 
163 (2009) 

[51] P. E. Theodorakis, W. Paul, and K. Binder, EPL 88, 
63002 (2009) 

[52] H.-P. Hsu, K. Binder, and W. Paul, Phys. Rev. Lett. 103, 
198301 (2009) 

[53] H.-P. Hsu, W. Paul, S. Rathgeber, and K. Binder, Macro- 
molecules 43 1592 (2010) 
[54] H.-P. Hsu, W. Paul, and K. Binder, Macromolecules 43, 

3094 (2010) 

[55] P. E. Theodorakis, W. Paul, and K. Binder, Macro- 
molecules 43, 5137 (2010) 

[56] P. E. Theodorakis, W. Paul, and K. Binder, J. Chem. 
Phys. 133, 104901 (2010) 

[57] Bates F.S., Fredrickson G.H., Annu. Rev. Phys. Chem. 
41, 525 (1990). 

[58] I.Ya. Erukhimovich, A.R. Khokhlov, Polymer Sci. 35, 
1522 (1993). 

[59] K. Binder, Advances in Polymer Science, 112, 181 
(1994). 

[60] I. Erukhimovich, Chapter 11, in: Nanostructured Soft 
Matter: Experiment, Theory, Simulation and Perspec- 
tives, NanoScience and Technology, edited by A. V. 
Zvelindovsky, Springer, New York, 2007. 

[61] I. W. Hamley, The Physics of Block Copolymers (Oxford 
Univ. Press., New York, 1998) 

[62] G. J. A. Sevink, A. V. Zvelindovsky, J. G. E. M. Fraaije, 
and P. Huinink, J. Chem. Phys. 115, 8226 (2001) 

[63] X. He, M. Song, H. Liang and C. Pan, J. Chem. Phys. 



26 



114, 10510 (2001) 
[64] Y. Wu, G. Cheng, K. Katsov, S. W. Sides, J. Wang, J. 

Tang, G. H. Fredrickson, M. Moscovits, and G. D. Stucky, 

Nat. Mater 3, 816 (2004) 
[65] K. Shin, H. Xiang, S. I. Moon, T. Kim, T. J. McCarthy 

and T. P. RusseU, Science 306, 76(2004) 
[66] H. Xiang, K. Shin, T. Kim, S. I. Moon, T. J. Mc- 
Carthy, and T. P. Russell, Macromolecules 37, 5660 

(2004); ibid38, 1005 (2005) 
[67] Y. Sun, M. Steinhart, D. Zschech, A. Rameshwar, G. H. 

Michler and V. Gosecke, Macromol. Rapid Commun. 28, 

369 (2005) 

[68] W. Li, R. A. Wickham, and R. A. Garbary, Macro- 
molecules 39, 806 (2006) 

[69] J. Feng and E. Ruckenstein, Macromolecules 39, 4899 
(2006) 

[70] J. Feng and E. Ruckenstein, J. Chem. Phys. 125, 164911 
(2006) 

[71] B. Miao, D. Yan, C. C. Han, and A.-C. Shi, J. Chem. 

Phys. 124, 144902 (2006) 
[72] B. Yu, P. Sun, T. Chen, Q. Jin, D. Ding, B. Li, and A.-C. 

Shi, Phys. Rev. Lett. 96, 138306 (2006) 
[73] P. Chen. X. Hi, and H. Liang, J. Chem. Phys. 124, 

104906 (2006) 
[74] Q. Wang, J. Chem. Phys. 126, 024903 (2007) 
[75] I. Erukhimovich and A. Johner, EPL 79, 56004 (2007) 
[76] G. J. A. Sevick and A. V. Zvelindovsky, J. Chem. Phys. 

128, 084901 (2008) 
[77] L. Leibler, Macromolecules 13, 1602 (1980) 
[78] I. Ya. Erukhimovich, Polymer Sci. USSR 24, 2223; 2232 



(1982) 

[79] LM. Lifsitz, A.Yu. Grosberg and A.R. Khokhlov, Rev. 
Mod. Phys. 50, 683 (1978); A. Yu. Grosberg and A. R. 
Khokhlov, Statistical physics of macromolecules, Ameri- 
can Institute of Physics, New York, 1994. 

[80] J.H. Hildebrand, R.L. Scott, The Solubility of Nonelec- 
trolytes, 3rd ed.\ American Chemical Society Monograph 
Series; American Chemical Society: Washington, DC, 
1950. 

[81] LYa. Erukhimovich, V. Abetz, R. Stadler, Macro- 
molecules 30, 7435 (1997) 

[82] P.G. de Gennes, Scaling Concepts in Polymer Physics 
(Cornell University Press, Ithaca, 1978). 

[83] LYa. Yerukhimovich, Polymer Sci. U.S.S.R. 21 470 
(1979); 

[84] LYa. Yerukhimovich, Letichii B.A. Polymer Sci. U.S.S.R. 
21 1394 (1979). 

[85] J.S. Higgins, H.C. Benoit, Polymers and Neutron Scat- 
tering. Oxford: Clarendon Press, 1996. 

[86] A.V. Dobrynin, LYa. Yerukhimovich, Macromolecules, 
26, 276 (1993), supplemented material available. 

[87] |http://www.gromacs.org| 

[88] W. F. van Gunsteren and H. J. Berendsen, Mol. Simul. 
1, 173 (1988) 

[89] G. S. Grest and M. Murat, Macromolecules 26, 3108 
(1993) 

[90] A. B. Bhatia and D. E. Thornton, Phys. Rev. B 2, 3004 
(1970) 

[91] F. H. StiUinger, J.Chem.Phys. 38, 1486 (1963) 



